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Abstract. A n X n matrix A has normal defect one if it is not normal, how- 
ever can be embedded as a north-western block into a normal matrix of size 
(n + 1) X (n + 1). The latter is called a minimal normal completion of A. A 
construction of all matrices with normal defect one is given. Also, a simple 
procedure is presented which allows one to check whether a given matrix has 
normal defect one, and if this is the case — to construct all its minimal nor- 
mal completions. A characterization of the generic case for each n under the 
assumption rank(A*j4 — AA*) = 2 (which is necessary for A to have normal 
defect one) is obtained. Both the complex and the real cases are considered. It 
is pointed out how these results can be used to solve the minimal commuting 
completion problem in the classes of pairs ofnxn Hermitian (resp., symmet- 
ric, or symmetric/antisymmetric) matrices when the completed matrices are 
sought of size (n + 1) X (n + 1). An application to the 2 X n separability problem 
in quantum computing is described. 



A matrix N E C nxn is called normal if N*N = NN*. For a non- normal 
A 6 C" xn it is natural to inquire what is the minimal p £ N such that A has a 
normal completion of size (n + p) x (n +p). In other words, what is the smallest p 
for which there exists a matrix of the form 



which is in fact normal? This minimal p is called the normal defect of A, denoted 
nd(yl), and a normal completion of size (n + nd(A)) x (n + nd(A)) is called minimal. 

The normal completion problem as above was introduced in [llj . and some ob- 
servations were made there. One of them is that among completions there 
exist those being scalar multiples of a unitary matrix. The smallest value of p 
required for such a completion is called the unitary defect of A, denoted ud(A), 
and the corresponding completions are called minimal unitary completions of A. 
Obviously, nd(A) < -ad(A). 

The latter inequality gives a simple upper bound for nd(A). Indeed, ud(A) is 
simply the number (counting the multiplicities) of the singular values of A different 
from \\A\\, and is therefore strictly less than n. Moreover, it was shown in [11] 
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that a minimal unitary completion of A can be constructed using the singular value 
decomposition (SVD) of A (see, e.g., [7] for the definition of the SVD). 

It is easy to find examples of matrices with nd(^4) different from ud(A). For 
instance, if A is normal and not a multiple of a unitary matrix then nd(A) = 
< ud(.A). However, in all such examples known until recently, the matrix A 
was unitarily reducible, that is, unitarily similar to a block diagonal matrix with 
non-trivial blocks of smaller size. It was therefore natural to ask pj] (see also a 
further discussion in [B]) whether the equality nd(^4) = ud(^4) holds for all unitarily 
irreducible matrices A G C nxn . We will show in Examples 12.231 and 12.241 that 
this question has a negative answer, and thus normal defect is indeed a different 
characteristic of a matrix then its unitary defect. 

A lower bound for nd(A) has been found in [8j: 

(1.2) nd(A) > m&x{i + {A*A - AA*),i-(A*A - AA*)}, 

where i+(X) and i-{X) denote the number of positive (negative) eigenvalues of a 
Hermitian matrix X. 

For 2x2 matrices, the unitary defect is at most 1. Therefore, any non- normal 
matrix A of size 2x2 has normal defect 1. Since trace(A*^4 — AA*) — 0, the 
righthand side of (|I.2| is if A G C 2x2 is normal and 1 otherwise. In other words, 
for n = 2 (jl.2p turns into an equality. We will show in Corollary 12.111 that this is 
also true for matrices of size 3x3. On the other hand, our Example 12. 121 reveals 
that for larger matrices a strict inequality in (|1.2[) is possible. 

The bulk of this paper is devoted to matrices with normal defect one. According 
to (|1.2p . all such matrices A must satisfy 

(1.3) rank(A*A - AA*) = 2 

which throughout the paper will be referred to as the rank condition. The manifold 
of n x n matrices satisfying (|1.3|) will be denoted 9Jt„. 

We obtain several equivalent characterizations of matrices A € C nx " with nd(A) 
1 . Among others, we prove in Section[2]that the following statements are equivalent: 

(i) nd(A) = 1. 

(ii) There exist a contraction matrix C € with ud(C) = 1, a diagonal 
matrix D G <C nxn , and a scalar ji € C such that 

A = CDC* + fil n . 

(iii) The rank condition holds, and the equation 

PA*(x\Ui + X2U2) — PA(x 2 ui +siu 2 ) 
has a solution pair X\,X2 G C satisfying 

N 2 - M 2 = d. 

Here u\, 1*2 are the unit eigenvectors of the matrix A* A — AA* corresponding to its 
nonzero eigenvalues \\ — d(> 0) and A2 = —d, and P = I n — U\U\ — U2U*, is the 
orthogonal projection of C™ onto rmll(A* A — AA*). 

(iv) There exist linearly independent x, y G C™ such that 

A* A - AA* = xx* - yy* 
and the vectors x,y,A*x,Ay are linearly dependent. 
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The equivalences (i)<S=>(ii), (i)-4=>(iii), and (i)<S==>(iv) are proved in Theorems 
12.1112. 31 and !2.4[ respectively. Note that version (ii) is better suited for construction 
of all matrices A with nd(A) = 1 (Section 12. ip . while (iii) is used to describe a 
procedure which allows one to check whether nd(A) = 1, and if this is the case — to 
construct all minimal normal completions of A (Section 12. 2[) . Finally, (iv) becomes 
handy when solving a separability problem originated in quantum computing (see 
description of Section [5] below). 

Section 12.31 provides a finer analysis which involves a certain basis construction 
and a lemma on symmetric and unitary extensions of certain matrices. This analysis 
allows us to refine the procedure from Section 12.21 and to describe the generic 
situation for each n under the assumption that rank condition (|1.3|) holds. In other 
words, certain topological characterization of the set of matrices with normal defect 
one in each matrix dimension is presented. 

We also obtain the following result: 

If A £ W ixn and nd(A) = 1 then one can choose a minimal normal completion 
of A to be real as well. 

It is natural to study the minimal normal completion problem also in the setting 
of real matrices, and we treat it in a separate Section [3l The real counterpart of 
the normal defect of a matrix A £ R nx ™, denoted rnd(A), is defined. The result 
stated above can be reformulated as follows: rnd(A) = 1 if and only if nd(A) = 1. 
(The question on whether rnd(A) = nd(A) for an arbitrary A £ M™ x ™ remains 
open.) This and other results in the real case are not immediate consequences of 
their complex counterparts, and required an additional study. Some statements in 
the real case are similar to their complex analogues, however there are also some 
differences. The real counterpart of characterization (ii) of matrices with normal 
defect one (see above) is obtained for matrices A £ R™ xn of even size n only, while 
such a characterization and a construction of all real matrices with rnd(A) = 1 in 
the case of odd n are left as an open problem. The real analogue of characterization 
(iii), as well as the procedure for verification that rnd(A) — 1 and for construction 
of all minimal real normal completions, have a slightly different form which splits 
into two cases. The generic situation in each matrix dimension is also described, 
however in the real case the analysis happens to be more straightforward than its 
counterpart in the complex case. 

In SectionHJ we show how to restate our results from Sections [2] and [3] in terms of 
commuting completions of a pair of Hermitian (resp., symmetric and antisymmet- 
ric) matrices, where the completed matrices are also Hermitian (resp., symmetric 
and antisymmetric). The results for pairs of Hermitian matrices are used then to 
solve an analogous problem in the class of pairs of symmetric matrices. 

In Section [5j we use the connection between the normal completion problem and 
the 2 x n separability problem, that was established in [12J . to obtain Theorem l5.ll 
which gives the easily verifiable necessary and sufficient conditions for a positive 
scmidefinite matrix M £ c2nx2n w jth a ran k one Schur complement to be 2 x n 
separable. Moreover, a new proof is given for the result by Woronowicz [13] (see 
Theorem 15. 2|) . which establishes, for n < 3, the 2 x n separability for a positive 
semidefinite matrix M £ C 2nx2n satisfying the Peres test. 

2. The complex case 
2.1. Construction of matrices with normal defect one. 
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Theorem 2.1. Let 4eC" x ™ be not normal. The following statements are equiv- 
alent: 

(i) nd(A) = 1. 

(ii) There exist a contraction matrix C G C txn with ud(C) = 1, a diagonal 
matrix D G C™ xn ; and a scalar /iGC such that 

(2.1) A = CDC* + fil n ■ 

(iii) There exist a unitary matrix V G C" x ™, a normal matrix N G C rix ™, and 
scalars t: <t < 1. /iGC such that 

(2.2) V*AV = MNM + fj,I n , 
where M = diag(l, . . . , 1, t). 

'A 



Proof. (i)<==>(ii) Let nd{A) = 1, and let 



G C(" +1 ) X (" +1 ) be a minimal 



normal completion of A. Then there exist a diagonal matrix A G C™ xn , a scalar 



/iGC, and a unitary matrix U 



U 2 i 



C/12 
U22 



G C( n+1 ) x (™ +1 ' such that 



'A 


X 






C/12" 




"A 


0" 




"TT* 


U*21 


y* 


z 




U21 


U 2 2_ 







/'_ 




TT* 
yi2 


C/22. 



(2.3) 

The latter equality is equivalent to the following system: 

(2.4) A = U 11 AU n +nU 1 2U^ = Un{A- f iI n )U n + f iI n , 

(2.5) x = U 11 hU^+ l xU 12 U22 = Ux 1 {K-^I n )U^, 

(2.6) y* = U 2 iAU n +fiU 2 2U^ = U 21 (A- nI n )U n , 

(2.7) x = C/ 21 AC7 2 * 1 +/xC/ 22 C722 = C/ 21 (A- M /„)C/| 1 + At . 

Setting C = Un and L> = A - fil n , we obtain |2~T|) from (f2~4|) . 

Conversely, if (|2.ip holds, we set Uu = C, A = D + fj,I n and obtain (I2.4[) . For 

tfn C/12 



C7 



C/21 U22 
z G C by (EU-llMl). Then 



a minimal unitary completion of C, we define x, y G C™ and 

.4 



holds, i.e., the matrix 



y 



£ C("+1)X(«+1) jt. 



a normal completion of A, and thus nd A = 1. 

(ii)^(iu) If (ii) holds, let C = F diag(l, ...,l,i)W* be the SVD of C (here 
V, W G C nx ™ are unitary, < t < 1, and M = diag(l, . . . , l,t) G C" XTl ). Then, 
clearly, = W*DW is normal, and (|2.2p follows. 

Conversely, if (|2.2p holds, then A^ = W*DW with D diagonal and iy unitary. 
Clearly, for C = V diag(l, . . . , 1, t)W* we have ud(C) = 1, and ([23]) follows. □ 



Remark 2.2. Observe that the matrix A given by 
and only if the product MNM is normal, that is 

(2.8) MNM 2 N*M = MN*M 2 Nt 
Since A^ itself is normal, (|2.8[) holds if and only if 

(2.9) MNZN*M = MN* ZNM. 
where Z = diag(0, ...,0,1). Partitioning A^ as 



happens to be normal if 



A 



A^o 9 
h* a 
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where a is scalar, and rewriting (|2.9|) block-wise, we see that it is equivalent to 

gg* = hh* , tah = tag. 

These conditions mean simply that g differs from h by a scalar multiple of absolute 
value one and, if to ^ 0, this scalar must be a/a. Consequently, A is not normal 
if and only if this is not the case. 

Observe also that if t ^ (so that M is invertible) and N is also invertible, then 
p.8p can be written as 

(2.10) M 2 N*N~ 1 = N- X N*M 2 . 

But N is normal, so that N* commutes with iV -1 . Condition (|2.10|) therefore 
means simply that N*N~ X (= N^ X N*) commutes with M 2 . In other words, A in 
this case is normal if and only if e„ := col(0, . . . , 0, 1) is an eigenvector of N*N~ 1 . 
In turn, this happens if and only if e n belongs to the sum of eigenspaces of N with 
the corresponding eigenvalues lying on the same line through the origin. 

Representation (|2.ip or (|2.2p in Theorem 12. 1[ together with Remark 12.21 allow 
one to construct all matrices A with nd(A) = 1. However, as we mentioned in 
Section [T] this does not give an easy way to check whether a given matrix has 
normal defect one. A procedure for this is our further goal. 

2.2. Identification of matrices with nd(A) = 1 and construction of all min- 
imal normal completions of A. In the following two theorems, we establish 
necessary and sufficient conditions for a matrix A to have normal defect one, and 
for any matrix A with nd(-A) = 1 we describe all its minimal normal completions. 
Here and throughout the rest of the paper, we set T = {z 6 C: \z\ = 1}. 

Theorem 2.3. Let A e C nxn . Then 

(i) nd(A) = 1 if and only if T&nk(A* A — AA*) = 2 and the equation 

(2.11) PA*{xiui +x 2 u 2 ) = PA(x 2 ui+xiu 2 ) 
has a solution pair Xx,x% £ C satisfying 

(2.12) \ Xl \ 2 - \x 2 \ 2 = d. 

Here U\, u% € C" are the unit eigenvectors of the matrix A* A — AA* corresponding 
to its nonzero eigenvalues Ai = d(> 0) and X 2 — —d, and 

(2.13) P = I n — U\U* — u 2 u 2 

is the orthogonal projection of C n onto null(^4*v4 — AA*). 

(ii) If nd(A) — 1, X\ and x 2 satisfy (|2.1ip and (|2 . 1 2[) . and [i G T is arbitrary 
then the matrix 

A li(x\U\ + x 2 u 2 ) 

~p(x 2 ul + XlU 2 ) z 

is a minimal normal completion of A. Here 

(2.15) z = an — — (x 2 (ai 2 xi — a 2 ix 2 ) + x\(ai 2 xi — a 2 \x 2 )) 
and 

(2.16) an = u*Aui, Qi2 = u*Au 2 , a 2 i — u 2 Aui. 
All minimal normal completions of A arise in this way. 



(2.14) B = 
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Theorem 2.4. Let A G C nxn . Then nd(^4) = 1 if and only if there exist linearly 
independent x,y G C" such that 

(2.17) A* A - AA* = xx* - yy* 

and the vectors x,y, A*x, Ay are linearly dependent. In this case, there exist z G C 
and v G T such that the matrix 

A vx 

y* z 



B = 



(2.18) 
is normal. 

In order to prove Theorems 12 .31 and 12 .41 we will need several auxiliary statements. 

Lemma 2.5. Let A G C nx ". Then xid(A) = 1 if and only if there exist linearly 
independent vectors i,i/6 C™ and a scalar z G C such that 

(2.19) A* A - AA* = xx* - yy*, 

(2.20) (A- zl n )*x = (A- zl)y. 

'A 



Proof. If nd(yl) = 1 then there exists a normal matrix B = 



y 



G c(™+ 1 ) x ( n + 1 ) 



The identity B*B = BB* is equivalent to (|2~T9| - (|2~20"1) (the identity x*x = y*y fol- 
lows from (|2.19p since tr&ce(A* A — AA*) = 0, and is therefore redundant). Clearly, 
x and y are linearly independent, otherwise the right-hand side of (|2.19p is and 
A is normal. 

Conversely, if x, y G C" are linearly independent, z G C, and (|2.19p - (|2.20[) hold 

then the matrix B — G c( n+1 ) x (' i+1 ' ) is normal. Since the right-hand side 

y z 

of (|2.19j) is not 0, the matrix A is not normal, thus nd(^4) = 1. □ 

Corollary 2.6. If nd{A) = 1 then rank(AM - AA*) = 2. 

The rank condition above is easy to check. Its equivalent form is that the char- 
acteristic polynomial of A* A — A A* can be written as 

(2.21) det(^*A - AA* - XI n ) = (-1)"A"- 2 (A 2 - d 2 ), 

with some d > 0. If this condition is satisfied, one can find the unit eigenvectors 
u\ and ui of the matrix A* A — AA* corresponding to its eigenvalues Ai = d and 
A2 = — d, which are determined uniquely up to a scalar factor. There is more 
freedom in a choice of other eigenvectors 1*3, . . . , u n , which form an orthonormal 
basis of null(yl*^4 — AA*). Suppose that such vectors are chosen. Then U — 



[ui 1t 2 1t 3 

satisfies 
(2.22) 
where 
(2.23) 



G C nx ™ is a unitary matrix, and the matrix A — U*AU 
A* A - A A* = H, 

i? = diag(d,-d,0,...,0)GC" x ™. 



Lemma 2.7. If A £ C" xn satisfies (|2~22f then A has the form 
(2.24) A 



an ai2 

0-21 022 

w* q* 
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where aij = 1,2) are scalars, 

(2.25) an = a 2 2, 

and u*,v*,w*,q* E C n ~ 2 satisfy 

( c\ c\t~>\ * * * * * 

(2. 2oj uu = qq , vv = ww , uv = wq , uw = vq . 

Proof. We have 

trace(Xff) = tr&ce(AA*A - A 2 A*) = 0, 
which implies (]2.25[) . Then, from the equality 

(A* A- AA*) 12 =0 
we obtain that uv* — wq*. Next, from the observation 
tr&ce{A*HA + AHA*) = tmce(A*(A*A-AA*)A + A(A*A-AA*)A*) 

= trace(l* 2 2 2 - A*AA*A + AA*AA* - A 2 A* 2 ) = 

we obtain 

(2.27) uu* - vv* + ww* - qq* = 0, 
and from the equality 

trace((A*A- AA*) 33 ) = 

we obtain 

(2.28) uu* + vv* = ww* + qq* . 

Combining (|2.27p and (|2.28p . we obtain that uu* = qq* and vv* = ww* . From the 
observation 

tT&ce(AHA) = trace(2(!*2-X4*)l) 

= trace(i!*2 2 - A 2 A* A) = 
we obtain that uw* — vq*. □ 
Lemma 2.8. Suppose that 



(2.29) 



A = 



an ai2 u 
a 2 i an v 
w* q* S 



A x 

y* z 



satisfies (|2.22p with H given by (|2.23p . Then the matrix B = 
is normal if and only if 

(2.30) a; = ool(xi,a;a,0,... ) 0) € C n , y = col(y u y 2 , 0, . . . ,0) e C n , 

(2.31) \xi\ 2 - \x 2 \ 2 = d, 

(2.32) yi=e te x 2 , y 2 = e i0 x u 
for some 9 € R, and the following identities hold: 

(2.33) (an - z)xi + a 2 ix 2 = (an - z)yi + ai 2 y 2 , 

(2.34) ai 2 xi + (an - z)x 2 = a 2 iyi + (an - z)y 2 , 

(2.35) u*xi+v*x 2 = w*yi+q*y 2 . 



g (J^(n+l)x(ri+l) 
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Proof. It follows from Lemma [23] applied to the matrix A as above (see also Lemma 



I2.7l which justifies that 022 = an) that nd(^4) = 1 and B 



g (£(n+l)x(rH-l) 



A x 

\y* z . 

is a minimal normal completion of A if and only if x and y are linearly independent 
and <|2.19|1 - (|2.20|) hold with A replaced by A. Since in this case H = xx* — yy*, 
the vectors x and y have the form (|2.30| . Indeed, for any vector h £ C™ which is 
orthogonal to y and not orthogonal to x, we have 

^ Hh = xx* h £ range(ff) n span(x). 

Similarly, for any vector g £ C n which is orthogonal to x and not orthogonal to y, 
we have 

7^ Hg = —yy*g £ range(ff) n span(y), 

thus both x and y are in range(-ff). Next, the identity H = xx* — yy* holds if and 
only if 

Nil 2 - |x 2 | 2 =d = \y 2 \ 2 - \yi\ 2 , x{x 2 = yiy 2 , 
or equivalently, (|2.32[) holds with some £ R. Clearly, (|2.20|) with A replaced by 



A, is equivalent to (|2~35|) - (|2~35|l . 

Remark 2.9. If A is as in Lemma 12.81 and B = 
minimal normal completion of A then so is 



A x 

y* z 



g (£(n+l)x(n+l) 



□ 



is a 



" A 


[IX 




In 


0" 




A x 




In 


0" 


py* 


z 











y* z_ 







/'. 



for any \i £ T. Therefore, if x, y and z are as in Lemma 12.81 then the matrix 



an 


012 


u 


a?\ 


an 


V 


w* 


<f 


S 


e~ l8 / 2 x 2 


e~ w ' 2 x l 






e -»e/2 
p -*e/2 





is a minimal normal completion of A. This observation leads to the following 
statement. 

Lemma 2.10. Suppose that 

an ai2 u 
A = a 2 i an w 
w* q* S 



satisfies (|2~22|) with H given by ([2~23]) . TTiera nd(A) 
Xi,x 2 € C satisfying (|2.12[) and smc/i i/iai 

(2.36) 

In £/iis case, i/ie matrix 



1 i/ and on/?/ if i/iere exist 



it Xi + v*x 2 = w*x 2 + q*x\. 



(2.37) 



an ai2 u xi 

a 2 i an v x 2 

w* q* S 

X2 Xi Z 



whe 



given by (|2 . 15[) . is a minimal normal completion of A. 
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Proof. The statement follows from Lemma 12. 8[ Remark 12.91 and the observation 
that, if yi = x 2 and y 2 = Wi, then (j2"lS5"j) becomes (|2~3g]) . Solving (|235|) or (pIM)) 
(which are equivalent in this case) for z gives (12.151) . □ 



Proof of TheoremUnJl (i) By Corollary[2~Hlthe rank condition, rank(A* A- AA*) = 
2, or equivalently, (|2.21[) with some d > 0, is necessary for A to have normal defect 
one, thus we can assume that this condition holds. Let u\ and u 2 be the unit 
eigenvectors of the matrix A* A — AA* corresponding to its eigenvalues Ai = d and 
A2 = — d, and let U3, . . . , u n be an arbitrary orthonormal basis of null(A*^4 — 
AA*). Define a unitary matrix U = \u\ ... u„] S C nx " and an isometry U' — 
[u 3 ■■■ u n ] £ C nx( - n - 2 l Then 

(2.38) U'U'* = P, 

and the matrix A = U*AU has the form (|2.29[) . where the scalars defined 
by (EH]), 

(2.39) u = u* x AU' 

According to Lemma \2. 10) nd(^4) 
is satisfied with some x\,X2 G 
be written as 



v = u*,AU' , 



w* = U'*Au\, q* = U'*Aii2- 



1 (and hence nd(A) = 1) if and only if ([2~36]) 
subject to (j2~T2"|) . By ([23^]) . equation (|23S|) can 



{U'*A*u l )xi + (U'*A*u 2 )x 2 = (U'*A Ul )x 2 + {U'*Au 2 )x 1 . 

Multiplying on the left by U' and taking into account that U' : C™ -2 
isometry satisfying (|2.38p . we obtain an equivalent equation 

(2.40) u*xi + v*x 2 = w*x 2 + q*X! , 
with the vectors 

(2.41) u*=PA* Ul , v*=PA*u 2 , w* = PA Ul , q* = PAu 2 

Note that these vectors are defined independently of the choice of u 3 , . 
(|2.40p is equivalent to (|2.1ip . this proves part (i) of this theorem. 

(ii) If nd(A) = 1 and A is defined as in part (i), then nd(A) = 1. 
[2~TU1 for any x 1: x 2 € C satisfying (f2~TT|) (or equivalently, pT36]) ) and (|2~T2]) . the 
matrix in (|2.37p is a minimal normal completion of A. By Remark 12.91 so is 



, u n . Since 
By Lemma 



(2.42) 



D 



an 
n 21 
w* 
~px 2 



ai2 
an 

q* 

~pxi 



u 

V 

S 




\1X\ 

HX 2 


z 



with an arbitrary /1 € T. By Lemma l2.8l all minimal normal completions of A arise 
in this way. Since A = U*AU and U is unitary, all minimal normal completions of 
A have the form 



B 



~u 


0" 


B 


~U* 


0" 







1 





1 





.4 



Ll(xiUi + X 2 U 2 ) 

z 



~p(x 2 ul + X\U 2 

where B is any of the matrices defined in (|2.42p . This proves part (ii) . 

Proof of Theorem \2.4\ If nd(A) = 1 then there exists a normal matrix 

'A x 

y* * 



□ 



B = 



£(n+l)x(n+l) 



10 D. S. KALIUZHNYI-VERBOVETSKYI, I. M. SPITKOVSKY, AND H. J. WOERDEMAN 



Ay — zy, i.e., x, y, A*x, Ay are linearly 



By Lemma [231 flUHl) holds and A*x-zx 
dependent. 

Conversely, suppose that A* A — AA* = xx* — yy* is satisfied with some linearly 
independent vectors x, y € C™, and the vectors x, y, A*x, Ay are linearly dependent. 
Clearly, in this case rank( J 4* J 4 — AA*) = 2. Choose an orthonormal eigenbasis 
U\, . . . ,u n of the matrix A* A — AA* and define a unitary matrix U = [ui . . . u n ] 
as in the proof of Theorem 12.31 By Lemma 12. 71 the matrix A = U*AU has the 
form (|2.29p . Define new linearly independent vectors x — U*x, y — U*y. Then 
A* A - A A* = xx* - yy*. As in the proof of Lemma [2T51 we conclude that (|2.30j) - 
(12.321) hold with x and y in the place of x and y. Since x, y, A*x, Ay are linearly 
dependent, so are x, y, A*x, Ay, i.e., the matrix 



x y A*x Ay 



x x y 1 auxi + a 21 x 2 ai 1 y 1 +a 12 y 2 
x 2 y 2 ai 2 xi+anx 2 a 2 ij7i + any 2 
u*x\ + v*x 2 w*yi + q*y 2 



has rank less than 4. Therefore u*x\ -\-v*x 2 and w*y\ +q*y 2 are linearly dependent. 
The identities yi = e l9 x 2 and y 2 = e %6 x\, together with the first three identities in 
(|2.26p , imply that there is G M such that 

u*xi + v*x 2 = {w*yi + q*y 2 )e 1 ' 1 ' = (w*x 2 + q*^^ 9 ^ . 

Putting x( = x\e l_ 2~ and x 2 = x 2 e~ l ^~ , we obtain 

* o i * o * — o i * o 
u x\ + v X2=w x 2 +q x\. 



By Lemma [2~TU1 nd(A) = 1, and therefore nd(A) = 1. 

The last statement of the theorem is obtained as follows. Let 



and 



col« x° 2 , 0, . . . , 0), y° = col(x°,3??, 0, . . . , 0), 



= an - - (x 2 {ai 2 ~x*{ - 021^2) + x\{ a i2x\ - 021^2)) 



(see (|2.15j) V By Lemma [2~T0l the matrix B = 
completion of A. Then 

B 



A 
,,0* 



X 

z 



is a minimal normal 



V 




' A x ' 




V* 




'A e-^x 


- <b-6 

e l T-_ 




0* 

V z _ 




■ e-d> 

e l ^ SL _ 




y* z 



□ 



is a minimal normal completion of A, i.e., we obtain (|2.18p with v = er 1 ^ . 

Applying Theorem 12.41 to 3x3 matrices, we obtain the following. 

Corollary 2.11. A matrix A £ C 3x3 has normal defect one if and only if rank( A* A 
AA*) = 2. Thus, for any 3x3 matrix A, one has 

nd(A) = max{i+(A*A - AA*),i_(A*A - AA*)}, 

i.e., (jl.2p in this case turns into the equality. 



Proof. The necessity of the rank condition has been established in Corollary 12.61 
The sufficiency follows from Theorem 12.41 since any four vectors in C 3 are linearly 
dependent. The second statement is obvious in the case where A* A — AA* = 0. It 
easily follows from the first statement in the case when r&nk(A*A — A A*) = 2. If 
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r&nk(A*A - AA*) = 3 then nd(A) = ud(A) = 2, and since A* A - AA* has three 
nonzero eigenvalues (counting multiplicities), we also have 

max{i+{A*A - AA*),i-{A* A - AA*)} = 2. 

This covers all the possibilities, since rank(A*A — AA*) ^ 1 due to the identity 
tvace(A* A - AA*) = 0. □ 



In the following example, we show that for n > 3 the rank condition (|1.3|) is not 
sufficient for nd(A) = 1. 

Example 2.12. Let 



A = 



Note that A (= A) is already of the form (T2T2T))) . We have 

"2 0" 

A* A - AA* = 



"0 





1 


—i 


2 















1 

—i 


1 

s/2 
I 

V2 


i 



-2 











Equation (|2.36p in this case takes the form 



Xi 



1" 




' 1 " 




= Xi 




i 


—i 



and it has no solutions with 



Xi\ 



X2 



= 2 > 0. Thus, by LemmaHini nd(A) > 1. 



Remark 2.13. If rank(A*yl — AA*) = 2 and tii, U2 £ C n are the unit eigenvectors 
of A*A — AA* corresponding to the eigenvalues Ai = d(> 0) and A2 = — d, then the 
vectors x = \fdu\ and y = \[du% satisfy ((2T7j). Indeed, u\ and 112 are orthogonal, 
hence linearly independent, and span(wi,M2) = range(A* A — AA*). For arbitrary 
a, b £ C, we have 

(A* A — AA*)(aui + bu2) = d{au\ — bu2) — d(u\u\ — U2U2){au\ + &M2), 

therefore A* A — AA* = d{u\u\ — u 2 u*,). 

However, as the following example shows, these x and y do not necessarily satisfy 
the conditions of Theorem [ 

Example 2.14. Let 



A = 



"0 





1 


i 












1 


i 


1 


1 






V2 


2 


2 ' 


i 


i 




s/3 




V2 


2 1 


2 



Then 



A* A ~ AA* = 



e\e x - e 2 e 2 , 



10 
0-100 



where e% and e 2 are standard basis vectors, which are the eigenvectors of the matrix 
A* A — AA* corresponding to its eigenvalues Ai = 1 and A2 = —1. However, the 
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arc 



vectors x = ei, y = e2, ^4*; 
linearly independent. 

On the other hand, we have nd(A) = ud(A) = 1: one of minimal normal com- 
pletions of A (in fact, its minimal unitary completion) is 



B 



We are now able to describe a procedure to determine whether nd(A) = 1 for 
a given matrix A £ C nxn , i.e., whether equation (|2.11[) in Theorem 12.31 has a 
solution pair xi,X2 £ C satisfying (|2 . 12|) . Moreover, this procedure allows to find 
all such solutions, and then, applying part (ii) of Theorem 12. 3i all minimal normal 
completions of an arbitrary matrix A with nd(-A) = 1. 



" 





1 


i 


V2- 








1 


i 


-1 


1 


1 











2 


2 l 


i 


i 




V3 







V2 


2 ' 


2 


-1 


V2 








. 



The procedure 



Begin 



Step 1. Verify the condition rank(A*j4 - A A*) = 2, or equivalently, (|2.21|) . If 
it is satisfied - go to Step 2. Otherwise, stop: nd(-A) > 1. 

Step 2. Rewrite (|2.11j) in the form (|2.40|) . where u, v, w, q are defined in (|2.41|) 
(see Theorem 12.31 for the definition of P, u\, and u-z). Let 



u = ur + iui, v = v R + iv/, q = q R + iqj, w = w R 
where u^, uj, v R , vj , q^, qj , wj, ~wj £ R™, and let 

xi = xm+ixn, x 2 =x R2 + ix I2 , 
where x R i, xn, x R2 , xi 2 £ Then (|2.40p becomes 



(2.43) 



Denote 



Q = 



q£ 



-u? 



q£ 
fqj 



q£ 



x n 

XR2 

X12 



qj 
q£ 



Find m = rank(Q). 

Step 3. Depending on to, consider the following cases. 

(1) to = 0. In this case, u = v = q = w = 0, and then (|2.40[) holds with any 
xi,x 2 £ C such that |a;i| 2 — |x 2 | 2 = d. Therefore, nd(A) = 1. Go to Step 4. 

(2) 1 < to < 3. In this case, ([2"13j) has nontrivial solutions. Let F £ R 4x ( 4 -" 1 ) 
be a matrix whose columns are linearly independent solutions of (12.431) . 



Then x = col(x R i, xn, x R2 , x I2 ) £ 
if x = Fh, with h £ M 4 " m . Setting F = 

U'2 



4 is a solution of (|2.43[) if and only 
Fi" 
F 2 



where Fi,F 2 £ 



d2x (4— m) 



write the condition \xi\ 2 — \x2\ 2 > as 



h T {¥ T l ¥ 1 -¥l¥ 2 )h > 0. 
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Therefore, nd(A) = 1 if and only if the matrix K = Ff Fi — F^F2 has at 
least one positive eigenvalue. If this is not the case, stop. Otherwise, for 
any h in the level hyper-surface ft T K/i = d, define 



xri 
xn 



= F 1 h, 



XR2 
XI2 



= F 2 h, 

xr2 + ixi2 satisfying (|2.40[) and such 



and thus obtain x\ — xr\ +ixn, x 2 = 
that \x x \ 2 - \x 2 \ 2 = d. Go to Step 4. 
(3) m = 4. In this case, (|2.43|) . and hence (|2.40[) . has no nontrivial solutions, 
and nd(A) > 1. Stop. 

Step 4. For each pair xi,x 2 £ C obtained in Step 3, find minimal normal 
completions of A as described by (12.14[) - (|2.16[) . 

End 

Remark 2.15. If m = 1 then K always has a positive eigenvalue and nd(A) = 1. 
Indeed, in this case F is a full rank matrix of size 4x3. Since null(F2) ^ {0} and 
null(F) = {0}, for a nonzero vector h G null(F2) we have 



h T Kh : 



h T ¥l Fi/i > 0. 



Example 2.16. Consider the 4x4 shift matrix 

"0100 

A 



Being a single Jordan cell, A is unitarily irreducible. Since 

-1 0" 

A* A -AA* = 





1 



A is not normal, and the rank condition is satisfied. Clearly, nd(yl) = ud(yl) 
and a minimal unitary (and thus normal) completion of A is given by 



(2.44) 



with any p G T. It is less obvious that any minimal normal completion of A 
has this form, as shown in [H Proposition 2]. We will give here an alternative 
explanation. Observe that d = 1, and the nonzero eigenvalues of A*A — AA* are 
Ai = 1, Aa = —1. Let efc, k = 1, . . . ,4, be the standard basis vectors in C 4 . Then 
u\ = e4 and u 2 = e\ are the eigenvectors of A* A — AA* corresponding to Ai = 1 
and A2 = — 1. We have 



P = I4 — e^e\ — e\e'i = 
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Since 

u* = PA* Ul = 0, 
we have 



= PA*u 2 = e 2 , w* = PAu x = e 3 , q* = PAu 2 = 0, 



l R 



T 



T 



e 2) 



e3, 



and then 



Q 





















-1 


1 
1 




Since m = rank(Q) = 2, all the solutions x = co\(xri,xii,Xr2-,%R2) £ K of the 
equation Qa; = are given by x = Fh where 



Fi 
F 2 



1 

1 







and h € M 2 is arbitrary. The matrix K = FjFi — F^F 2 = I 2 is positive definite, 
and therefore nd(A) = f . Define, for an arbitrary /i£l 2 such that h\ 



hl = l, 



XRl 
XII 



F1/1 



XR2 
XI2 



F 2 h = 0, 



and we obtain x\ — Xm +ixn = hi +ih 2 and x 2 — xr 2 + ixi 2 = satisfying (|2.40[) 



and such that |a;i| 



\X2\ 



= 1. Then we calculate 



an = u\Au\ = U2AU2 = 0, a\ 2 = ulAu 2 = 0, 021 = u^Aui = 0, 

z = an - - {x 2 (ai 2 xi - a 2 ix 2 ) + xi(a 12 xi - a 2 ix 2 )) = 0. 
a 

Thus, all minimal normal completions of A have the form 



A \ix\u\ 
JIX1U2 



10 
10 
1 

0000c 

p 



where C = fJ-xi and p = ~px\, with some fi G T, i.e., £ and p are arbitrary complex 
numbers of modulus one. We conclude that all minimal normal completions of A 
have the form (|2.44[) . i.e., are minimal unitary completions of A. 

We note that a similar argument can be applied to weighted shift matrices of 
any size n > 4, with all the weights of modulus one, thus recovering the result of 
Proposition 2 in [8]. We also note that, as was mentioned in [8], for n = 2 or 3 there 
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exist non-unitary minimal normal completions of such weighted shift matrices. Our 
procedure gives the full description of these completions B. Namely, for n = 2 

1 



B = 









~fiX2 



[IX 2 
fJ,Xi 



X\X2 



where /ieT and x\,x% £ C: \x\\ 



B = 



| #2 1 2 = 1 are arbitrary, and for n = 3 



" 


1 





/J,X 2 








1 














HXi 


~pxi 





Jlx 2 






where /i £ T and x\ £ C, X2 £ K: |^i| 2 — x\ — \ are arbitrary. 

Later on, we will present more examples of application of this method, see Ex- 
amples (j2~2"3) and $TZty . 



2.3. The generic case. The procedure described in Section [2721 which is based 
on part (i) of Theorem 12. 3[ allows to check whether a given matrix A £ C" x " has 
normal defect one, and if this is the case — to solve the system of equations (|2.11| - 
(|2.12| . Part (ii) of Theorem 12.31 describes all the minimal normal completions of 
A. That procedure verifies first the necessary rank condition, and then uses only 
the two nonzero eigenvalues, Ai = d and A2 = — d, and the two corresponding unit 
eigenvectors, u\ and 112, of A* A — AA* . The vector equation (|2.43|) appearing in 
that procedure is equivalent to a system of 2n real scalar linear equations with 4 
unknowns. 

In this section, we show how the procedure in Section 12.21 can be refined by 
using a special choice of the eigenbasis for the matrix A*A — AA*, i.e., a special 
construction of orthonormal eigenvectors U3, . . . , u n corresponding to the zero 
eigenvalue. This additional analysis is rewarded by obtaining a system of n — 2 
(as opposed to 2n) real linear equations with 4 unknowns. Moreover, it allows us 
to describe the generic situation under the assumption that the rank condition is 
satisfied. The refined procedure is based on the following theorem (the proof of 
which is given later in this section). 

Theorem 2.17. Let A £ C nxn satisfy the rank condition (|1.3|) (or equivalently 
(|2.21l) ), and let u\ and U2 be the unit eigenvectors of the matrix A* A — AA* corre- 
sponding to its eigenvalues Ai = d(> 0) and A2 = —d. Then 

(i) There exist orthonormal vectors U3, ...,u n £ rmll(A* A — AA*) (and thus 

the matrix W — \u\ . . . u n ~\ £ C" x ™ is unitary) such that the matrix 

A = W*AW has the form 



(2.45) 



A = 



an a\2 u 
Q2i an v 

v 



T u T S 



with aij 's defined in (|2.16p . 
(ii) nd(A) — 1 if and only if the equation 

(2.46) Im(u*xi + v*x 2 ) = 
has a solution pair Xi,X2 £ C satisfying 

(2.47) \ Xl \ 2 - \x 2 \ 2 = d. 
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(iii) If nd(A) — I, x\ and xi satisfy p.46|) and ()2.47[) . and [i G T is arbitrary 
then the matrix B defined in (12.14|) is a minimal normal completion of A. 

All minimal normal completions of A arise in this way. 

Remark 2.18. The matrix W in Theorem 12.171 can be constructed explicitly as 
will become clear from the proof of the theorem. 

Let A G C nxn satisfy the rank condition. We define the vectors u*,v*,w*,q* G 
range(P) C C™ by (|2.41[) (see also Theorem [23] for the definition of P, u\, and 1*2). 
Since these vectors can be viewed as the images of vectors u*,v*,w*,q* G C"~ 2 
under an isometry so that (|2.26j) holds (see Lemma 12.71 and the proof of Theorem 
I2.3|) . we have 

(2.48) uu* = qq*, vv* = ww*, uv* = wq*, uw* = vq*. 
The first three equalities mean that the linear operator 

(2.49) X: span(u T ,v T ) — >span(q*,w*) 
defined via 

(2.50) X : u T 1 — > q*, v T 1 — ► w* 

is a well defined unitary operator. In order to interpret the last equality in (|2.48[) 
we need an intermission for some definitions and results on (complex) symmetric 
operators and matrices (see, e.g., Section 4.4]). 
For a subspace TL in C fc , denote its complex dual by 

TL:={heC k :heTL}. 

We will say that a C-linear operator L : TL — * TL is symmetric if h T Lg = g T Lh 
(or, equivalently, (Lg,h) = {Lh,g) in the standard inner product in C fc ) for all 
g,h G TL. It is clear that a C-linear operator L: C k — ► C fe is symmetric if and 
only if its matrix in a standard basis of C k is complex symmetric, i.e, L = L T . In 
general, a C-linear operator L: TL — > TL is symmetric if and only if its matrix in 
any pair of orthonormal bases B — {hj} k =1 C TL and B = {hj} k =1 C TL is complex 
symmetric, i.e., (Lhj, hi) = {Lhi, hj), i,j = l,...,k. 

We can restate this also in a coordinate-free form. Let Q and TL be two subspaces 
in C fe . For a C-linear operator L: Q — ► TL we define its transpose as the unique 
C-linear operator L T : TL — > Q which satisfies h T Lg = g T L T h (or, equivalently, 

(Lg,h) = (Lh,g)) for all g e G, h G TL). Explicitly, L T h = L*h for every h G TL. 
This, in particular, implies that (LM) T = M T L T for two C-linear operators L, M . 
If one interprets a vector h G TL as an operator h: C — > TL then its transpose h T can 
be interpreted as the operator h T : TL — > C, and then the identity h T Lg = g T L T h 
can be viewed also as g T L T h = (h T Lg) T . Finally, a C-linear operator L: TL — > TL 
is symmetric if and only if L = L T . 

We also observe that a matrix of a C-linear operator L: Q — > TL in a pair of 
orthonormal bases Si and $2 and a matrix of L T : TL — > Q in the pair of orthonormal 
bases £?2 and Si are transposes of each other. 

Lemma 2.19. Let TL and C be subspaces in C k such that TL C C, and letY : TL — > C 
be an isometry such that P^Y : TL —* TL is a symmetric operator. Then there exists 
a unitary and symmetric operator Y : C — * C such that Y\n = Y. 
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Proof. We have 

C = n®g®ic, 

where Q = range (P^q^F) and AC = C (TC © Q). The Takagi decomposition (see, 
e.g., [2]) of the symmetric operator PjjY, in a coordinate-free form, is 

PjjY = UY,U T 1 

where U : C dlm ^ — > H is a unitary operator such that E7ei, . . . , Ue dim ^ are the 
eigenvectors of PjjYP^Y, and E : C dlm ^ — » C dlm ^ w ^ is an operator whose matrix 
in the standard basis of C dlm ^^ is diagonal, with the singular values of PjtY on 
the diagonal. The operator PgY can be represented as 

PgY = V{l CiiMH) ~ Z 2 ) 1/2 U T , 

where V : C dlm ' w ^ — > Q is a coisometry, with 

^lra„ g c(/ cdlm(T1) -£2) : range(/ C d im( «) - E 2 ) -> Q 

unitary. Define the operator J : JC — > AC for some pair of orthonormal bases 23 = 
{^}f=i i>C) CK,B= {^}fZ {K) C AC as 

dim(/C) 

a J K j I = X! a J^'' a i> • • • > a dim(K) e c. 




3=1 



Clearly, J is symmetric, and the matrix of J in the pair of bases B and B is -Tdim(JC) ■ 
It is straightforward to verify that the operator 

H H 

XJYU T 17(/ C di«(«) - E 2 )V2yT 0] © © 

y := y(/ C din, ( «) - Y?fl 2 U T —VYiV T o ■. g ^ g 

jj © © 

AC AC 

has the desired properties. □ 

Corollary 2.20. The unitary operator X defined by (|2.49[) (|2.50|) can be extended 
to a unitary and symmetric operator X: rangc(P) — * range(P). 

Proof. Since u T ,v T S range(P) and q*,w* S range(P) (see (|2.4ip ). the operator 
X can be viewed as an isometry X : span(u T , v T ) — > range(P). The last identity 
in (|2.48jl means that the operator 

-P sp an(u-,v*)^: span(u T ,v T ) -> span(u*,v*) 

is symmetric. Then the statement of this corollary follows from Lemma[2J~9l where 
we set k = n — 2, 

Ti = span(u T ,v T ) = span(PA T ui, PA T v,2), 
Y = X and C = range(P). □ 

Remark 2.21. It can be shown that the unitary and symmetric operator X in 
Corollary 12.201 can be constructed bypassing Lemma 12.191 and using instead the 
following remarkable theorem from [13] : If A G C" x ™ and x, y S C" are such that 
A* A — AA* = xx* — yy* then there exists an antiunitary involution i on C n such 
that lx = y and lAl — A*. (A mapping i: C n — > C™ is called an antiunitary 
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involution if i 2 = I n and (ih, ig) = (g, h) for every h, g G C n , in the standard inner 
product in C n .) Our Lemma \2 . 1 91 seems to be of independent interest, and can be 
applied to other problems as well. 



Proof of Theorem \ 2. 17\ (i) The operator X in Corollary 12. 20( which is constructed 
as in Lemma 12.191 for the given matrix A, is symmetric and unitary, and thus has 
a Takagi factorization (see [7]) X = GG T where G: C™~ 2 — > range(P) is unitary. 
One can view G as an isometry G' : C™~ 2 — > C" with the same range as G: 

range(G') = range(G) = range(P). 

Clearly, the columns U3 , . . . , u n of the standard matrix of G' are orthonormal, and 
hence, together with u x and u 2 , form an orthonormal eigenbasis of A*A — AA*. We 
also have G'G'* = P. We then extend X to the operator X' = G'G' T : C n -> C™. 
The operator represented by the matrix A in the standard basis of C™, in the basis 
U\, . . . , u n has the matrix 



A = 

We have 
{ulAG'f 

and similarly, 



G" 



G 



an 

021 

G'*Aui 



ai2 
an 
G'*Au 2 



u\AG' ' 
u* 2 AG' 
G'* AG' 



= G 



IT 



A T ui = G'*G'G lT A T ui = G'*X'A T ui = G'*X'PA T ui 



= G'*PAu 2 = G'*Au 2 , 



u* 2 AG') T = G'*Aui. 



"3 



Setting u = u\AG' , v = u* 2 AG' , S = G'*AG', and W = [u x u 2 
we see that A = W*AW has the form (|2.45|) . which proves part (i). 

(ii) It follows from Lemma [2.101 that nd(A) = nd(A) = 1 if and only if there 
exist x x ,x 2 € C satisfying (|2.47|) such that 

u*xi + v*x 2 = v T x 2 + u T xi. 



□ 



Since the last equation is equivalent to (12.46|) . this proves part (ii). 
(iii) This part is proved in the same way as part (ii) of Theorem 12.31 

Let u, v G C™~ 2 be as in Theorem 12.171 u = ur + iui, v = Vr + where 
u^uJ^Vx^J G R"~ 2 . Let xi = x R i +ixn, x 2 = x R2 +ix I2 , where x RX , xn, x R2 , 
xj 2 G R. Then (|2.46p can be written as 

(2.51) Qx = 0, 

where 



(2.52) 



Q=[- 



>(n-2)x4 



XRI 

xn 

XR2 
Xll 



Remark 2.22. Observe that replacing u, v, w, and q in (|2.43[) by u, v, v, and u 
as in Theorem 12 . 1 71 we obtain an equivalent condition, i.e., equation (|2.51[) replaces 
(|2.43| with the matrix Q replacing Q. Thus instead of 2n real linear equations 
with 4 unknowns we obtain n — 2 real linear equations with 4 unknowns. Let 
m = rank(Q). Then, for all possible cases of m, the procedure for checking whether 
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nd(A) = 1, and if this is the case - for constructing all minimal normal completions 
of A, is exactly the same as described is Section |2~2"1 with Q in the place of Q. 

We describe now the generic situation for each n, under the assumption that 
rank(j4*yl — AA*) = 2. In other words, we obtain certain topological characteriza- 
tion of the set of matrices with normal defect one in each matrix dimension. 

The generic case 

Let A <G C" xn satisfy the rank condition. Consider the following possibilities for 
the value of n, and describe the situation for each case separately. 

n = 2 or n = 3: Vectors u, v as in Theorem 12.171 do not arise (when n = 2) 
or are scalars (when n = 3). Then m — rank(<5) < 1, and nd(A) = 1 (for 
the case where m = 1 it follows from Remark l2.15p . This gives a new proof 
of the statement on 2 x 2 matrices in Section [T] and of Corollary 12. Ill 

n = 4 or n = 5: In these cases, m < 2 (resp., m < 3). Thus, equation (|2.51j) 
(see also (|2.52p ) has nontrivial solutions. Both the situation where the 
matrix K, constructed from Q instead of Q, has at least one positive eigen- 
value (in which case nd(A) = 1) and where K has no positive eigenvalues 
(in which case nd(A) > 1) occur on sets with nonempty interior in the 
relative topology of the manifolds TI4 and WI5 (see page [2] for the definition 

ofart„). 

n > 6: In this case, generically m = 4, thus (|2.5ip has no nontrivial solutions. 
Therefore, generically nd(A) > 1. Still, there are matrices A with nd(A) = 
1, which can be constructed, e.g., using Theorem 12. II and Remark 12.21 



2.4. Normal defect and unitary defect. In this section, we give two examples 
which show that the question in [11] (see also [8]) asking whether nd(A) = ud(A) 
for any unitarily irreducible matrix A has a negative answer. In the first example, 
A has a single cell in its Jordan form, and in the second example, A has three 
distinct eigenvalues. We also present all minimal normal completions of A in both 
examples. 



Example 2.23. Let 



Then 



A = 





"2 





1" 




1 





f 


A*A = 





1 


1 


, AA* = 





2 


1 




1 


1 


2 




1 


1 


2 



and 



A* A- A A* 



1 

-1 




Clearly, the rank condition holds. It follows from Corollary 12.111 that nd(^4) = 1. 
Since the only eigenvalue of A is 1, and A — I is nilpotent of order 3, A has a single 
cell in its Jordan form, and hence A is unitarily irreducible. The characteristic 
polynomial of A* A is 

p(A) = (2-A) 2 (l-A) + 2A-3. 
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We have p(0) = 1 > 0, p(l) = -1 < 0, p{2) = 1 > 0, and p(4) = -7 < 0. Therefore, 
p(A) has three distinct roots, in intervals (0,1), (1,2), and (2,4), i.e., A has three 
distinct singular values. Therefore, ud(A) = 2. 

We also observe that A has the form (|2.45|) . The procedure described in Section 
together with Theorem 12.31 (or its refined version described in Remark 12.221 



together with Theorem 12. 17p gives that all minimal normal completions of A have 
the form 

" 1 

1 

1 

jlX2 JlXl 



B = 



fJ,Xi 

\xx 2 

1 



with arbitrary fi G T, and x\ G C, x 2 G 
Example 2.24. Let 



satisfying \x\ 



= 1. 



A 



1 

1 
1 



Changing the basis, we obtain A = U*AU, where 



U = 



is unitary and 



satisfies 



.4 



A* A - A A* 



- 





r 


1 

71 


1 

v/2 





I 

L V2 


t 

72 





' 3i 


i 


i ■ 


4 
7i 


4 
3/ 




4 

1 

. V2 


4 
1 

V2 






3 
0-3 




Clearly, mnk(A*A - AA*) = mnk(A*A- AA*) = 2. By Corollary[2T[J nd(A) = 1. 
We also observe that A has the form (|2.45|) . The matrix A is unitarily irreducible. 
Indeed, matrices 



Re(A) 



and Im(A) 



do not have common eigenvectors. 

Next we show that ud(^4) = 2. The characteristic polynomial of A* A, 



p(A) = (l-A)(2-A) 



13 



- A 



13 



A- 



17 

T' 



has values p(0) = § > 0, p(l) = -1 < 0, p(2) = | > 0, p(5) = -9 < 0. Therefore, 
p(A) has three distinct roots, in intervals (0,1), (1,2), and (2,5), i.e., A has three 
distinct singular values. Thus, ud(A) = 2. Note that in this example A has three 
distinct eigenvalues, Ai = i, A 2 = v ^f +i , A 3 = 
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The procedure described in Section 12.21 together with Theorem 12.31 (or its re- 
fined version described in Remark 12.221 together with Theorem I2.17j) gives that all 
minimal normal completions of A have the form 

H(hi + ih 3 ) 



B 



3i 

4 

— H 
4 

1 

V2 



4 
3i 
4 
1 
V2 



p(h 2 — iha) fi(hi + ihs) 



with arbitrary /j, G T, and hx^h^hs G K 
normal completions of A have the form 



B 



l 






M(^2 - ih 3 ) 


/ll/l3+5/l2?13+i(3-2/ll/l 2 +2/l2) 

3 ^~ 



"0 


1 





1 





1 





1 


2 l 





— hi+h 2 


— 2h 3 +i{h 2 -hi) 
V V2 



3. Correspondingly, all minimal 





, hi+h-2 



>" V2 

2h 3 +i(h 2 - 



* V2 
h 1 h 3 +5h 2 h3+i(3-2h 1 h 2 +2hl) 
3 



3. The real case 

nx ™. We define the real normal defect of A, rnd(A), as the minimal 

A 



Let A e 

possible nonnegative integer p such that a matrix 



l (n+p)x(n+p) jg normal 



(such a matrix with the minimal possible p is a minimal real normal completion of 
A). It is clear that rnd(A) > nd(^). We will show later (in Corollary [23) that 
nd(v4) = 1 if and only if rnd(A) = 1, while we leave open the question whether 
rnd(A) = nd(A) is always the case. 

We also define the orthogonal defect of A as the minimal nonnegative integer s 
A % 

such that a matrix ^ 6 jj( n + s )x(™+ s ) j s a multiple of an orthogonal matrix 

(such a matrix with the minimal possible s is a minimal orthogonal completion of 
A). In fact, the orthogonal defect of A coincides with ud(j4), so that it does not 
require a separate notation. Indeed, since a minimal orthogonal completion of a real 
matrix is obtained using the same construction as for a minimal unitary completion 
(see [H]), the only difference being that the real SVD is involved, the size of this 
minimal orthogonal completion is the same as for a minimal unitary completion. 

It is clear that rnd(A) < ud(A), and we will show later, in Example 13. 5[ that 
there are orthogonally irreducible matrices A such that the strict inequality takes 
place. 

3.1. Construction of real matrices of even size with real normal defect 

one. The following theorem is an analogue of Theorem [2J] for the case of real n x n 
matrices with n even. 



Theorem 3.1. Let A e R nxn , where n = 
statements are equivalent: 

(i) rnd(A) = 1. 

(ii) There exist a contraction matrix C G 
nal matrix D € R nx " of the form 



2k, be not normal. The following 
l nxn with ud(C) = 1, a block diago- 



(3.1) 



D = diag 



( 


ai 


lh~ 




at 


fit 






ai_ 


) * * • 3 




Oil 
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and a scalar \i £ K such that 

(3.2) A = CDC T + [il n . 

(iii) There exist an orthogonal matrix V €E R nx ™ 7 a normal matrix N £ 
and scalars t, fj, G R, wii/i < t < 1, suc/i £/ia£ 

(3.3) y T AV" = MNM + }il n , 
where M = diag(l, . . . , 1, t). 



Proof. (i)4=>(ii) Let rnd(yl) = 1, and let 



.4 



E (n+l)x(n+l) be a minima l 



real normal completion of A. Then (see, e.g., [H Section IX. 13]) there exist a block 
diagonal matrix A £ R™ x ™ of the form 







ft" 






Pi 


( 


-A 


Mi. 


, . . . j 







A = diag 
R, and an orthogonal matrix O = 



On 12 
O21 O22 



, A*2fc 



(n+l)x(n+l) 



a scalar fi G 
that 

(3.4) 



Here we used the fact that n is even, and thus the (n + 1) X (n + 1) real normal 
matrix T has at least one real eigenvalue. The last equality is equivalent to 
the following system: 



'A 


X 




'On 


12 ~ 




A 


0" 




\o u 




T 

V 


z 




O21 


2 2_ 











OI2 


022. 



(3.5) 


.4 


= O u AO£ 


+ ii0 12 T 12 = On (A - M-Wu + t*I n , 


(3.6) 


X 


= OnAO£ 


+ MO12O22 =O u (A-/x/ n )Oj 1 , 


(3.7) 


T 

y 


- 21 AO^ 


+ A t0 22 0f 2 = 21 (A-/z/ n )0£, 


(3.8) 


z 


= 21 AOj 1 


+ /i0^ 2 = 2 i (A - + /1. 


Setting 


C = O n 


and D = A - 


Ai/ Tl , we obtain (|3.2[) from (|3.5[) . 



o 



Conversely, if 
"On 12 



holds, we set On = C, A = D + and obtain 

Q mi 

2 i 2 2^ 

z G R by (JSS])-(1S3). Then (03) holds, i.e., the matrix 



5|). For 

a minimal orthogonal completion of C, we define x,y € R ra and 

(n+l)x(n+l) j g 



.4 



a real normal completion of A, and thus rnd(A) = 1. 

(ii)^(iii) If (ii) holds, let C = V r diag(l, . . . , l,t)W T be the SVD of C (here 
V,W e R" xn are orthogonal, < t < 1, and M = diag(l, . . . , 1, t) S R" xn ). 
Clearly, TV = W T DW is normal, and ([3?3l) follows. 

Conversely, if (|3.3[) holds, then N = W T DW with D block diagonal of the form 
(l3~Tj) and W orthogonal. For O = V diag(l, . . . , l,t)W T we have ud(O) = 1, and 
follows. □ 



Remark 3.2. Rcmark l2.2l can be restated in the real case as follows. The matrix A 
of even size, given by (|3.3|) . is not normal if and only if in the matrix N partitioned 

as 

'N g 
h T 



N 
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(where a is scalar) g ^ h and, in the case where to = 0, also g ^ —h- Moreover, if 
both M and N are invertible then A is not normal if and only if the standard basis 
vector e„ is not an eigenvector of N T N~ 1 . The statement in the last sentence of 
Remark 12.21 is. in general, not valid in the real case. 



Open problem: What is an analogue of Theorem 13. II for the case where n is 
odd? 



In the case of even n, similarly to the complex case, representation (13. 2|) or ()3.3[) 
in Theorem l3.11 along with Remark 13.21 allow one to construct all matrices A with 
rnd(A) = 1. However, this does not give a way to check whether a given real matrix 
has real normal defect one. A procedure for this is our further goal. 

3.2. Identification of matrices with rnd(A) = 1 and construction of their 
minimal real normal completions. In the following theorem, which is the real 
counterpart of Thcorcm l2.31 we establish the necessary and sufficient conditions for 
a matrix A s M™ x " to have real normal defect one, and for any A with rnd(A) = 1 
we describe all its minimal real normal completions. 

Theorem 3.3. Let A G R" x ". Then 

(i) rnd(A) = 1 if and only if r&nk(A T A — AA T ) = 2 and at least one of the 
equations 

(3.9) {PA T ui - PAu 2 )x l + {PA T u 2 - PAm)x 2 = 0, 

(3.10) (PA T ui + PAu 2 )x 1 + {PA T u 2 + PA Ul )x 2 = 0, 
has a solution pair xi,x 2 £K satisfying 

(3.11) x\-x\ = d. 

Here ui,u 2 € K™ are the unit eigenvectors of the matrix A T A — AA T corresponding 
to its nonzero eigenvalues Ai = d(> 0) and X 2 = —d, and 

(3.12) P = I n — uiu^ — u 2 u^ 

is the orthogonal projection ofW 1 onto null(^4 T v4 — AA T ). 

(ii) If rnd(A) = 1 then at least one of the following cases occurs: 

Case 1. If Xi and x 2 satisfy (|3.9[) and (|3.11[) then the matrix 



(3.13) B 1 



A xim + x 2 u 2 

x 2 uf + XiU 2 z 



is a minimal real normal completion of A. Here 

(3.14) z = an - -Axi + x 2 ){ai 2 xi ~ a 21 x 2 ) 

d 

and 

(3.15) an = ujAui, ai2 = ufAu 2 , a 2 \ — u 2 Aui; 

Case 2. If x\ and x 2 satisfy (|3.10[) and (|3.1ip then the matrix 



(3.16) B 2 



A xiu\ + x 2 u 2 

T T 
-x 2 u\ — X1U2 z 



is a minimal real normal completion of A. Here 

(3.17) z = a u + -(x 1 - aj 2 )(oi2a?i + 021X2) 

a 
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id (Xi 



defined by (|3.15[) . 



Any minimal real normal completion of A arises in this way, i.e., either as in 
Case 1 or as in Case 2 above. 



Proof. Since rnd(A) = 1 implies nd(A) = 1, it follows from Corollary 12.61 that 
the rank condition is necessary for rnd(A) = 1. For real A it takes the form 
rank(A T A — AA T ) = 2. Clearly, it is also equivalent to 



(3.18) 



det(A J A- AA 1 - AI„) = (-l)"A ,l " 2 (A 



d% 



with some d > 0. 

Without loss of generality, we can assume that the rank condition is satisfied. 
Then we find the unit eigenvectors u\,U2 G R" of the matrix A T A — AA T corre- 
sponding to its eigenvalues \\ = d and A2 = — d. Let U3, . . . , u n be an orthonormal 
basis for mi\\(A T A - AA T ). Then U' = [u 3 u n ] € R nx ( n ~ 2 ) is an isometry, 

and 



(3.19) 



U'u'T = p 



where P is defined in (|3 . 1 2|) . The matrix A = U T AU ', where U = \_u\ 
orthogonal, has the form 



is 



(3.20) 



A = 



an ai2 



an 

T 



(the identity an = ujAui = \1\Au2 is established in the same way as in Lemma 
2.7p . As in Lemma \2. 81 we obtain that rnd(A) = 1 (and therefore, rnd(A) = 1) if 
and only if there exist xi, x%, yi, j/2> z € K such that the matrix 



(3.21) 



B 



an 



021 

T 



an 

T 

q 

!)2 



u 

V 

S 




xi 

X2 



z 



»(n+l)x(n+l) 



is normal if and only if there exist x\, X2, yi, y-ii z € K such that 



(3.22) 
(3.23) 
(3.24) 
(3.25) 
(3.26) 



(an - z)xi + a 2 ix 2 
ai 2 xi + (an ~ z)x 2 

T T 
U Xi + V X2 

2 2 

x i -Vi 
XIX2 



(an ~ 
a 2 iyi 

T 

w yi 

vl 

2/12/2- 



z)yi 4 
f (an - 

l-<7 T 2/2, 
x\ = d, 



«122/2, 
" 2)2/2, 



It follows from (|3.25|) and (|3.26p that either yi = X2, 2/2 = xi or yi 
We will consider these two cases separately. 

Case 1: y\ = x 2 , 2/2 = Xi. Identities (|3.22() — ()3.25|) become 



-X2, 2/2 



-Xi. 



(3.27) 
(3.28) 
(3.29) 
(3.30) 



(an - z)xi + a 2 ix 2 = (an - z)x 2 + ai 2 xi, 
012X1 + (an - z)x 2 = a 2 ix 2 + (an - z)xi, 
(u T - q T )xi + (v T - w T )x 2 = 0, 



2 _ 2 _ 1 
x-\ x<y — a. 
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Clearly, (13.27|) and (|3.28l) are equivalent, and it follows from (|3.27p and (|3.30|) that 

a\2X\ - a 21 x 2 1 , . , , 

z = an = an — j(JCi + aigMaiaiBi - 021^2) 

xi — £ 2 a 



(cf. ([235]) ). Next, it follows from (pHO]) that 

m t = [/' t ^ t mi, w T = C/ /T A T u 2 , w T = t/' T A Ul , g T = [/ /T Au 2 . 

Multiplying both parts of these equalities on the left by U' and taking into account 
(13 . 19[) , we obtain vectors 

(3.31) u T = U'u T = PA T Ul , v T = U'v T = PA T u 2 , 

(3.32) w T = U'w T = PAui, q T = U'q T = PAu 2 . 

Since U' is an isometry, (|3.29|) is equivalent to 

(3.33) (u T - q r )xi + (v T - w T )x 2 = 0. 

Note that the definition of vectors u T , v T , w T , q T £ range(P) C R n in (|33Tj) - (j3~32"l) 
is independent of U', i.e., of the choice of basis vectors U3, . . . , u n in range(P) = 
nu \\{A T A-AA T ). 

Case 2: y\ = ~x 2 , y 2 = —x\. Identities (j3.22[) (j3.25|) become 



(3.34) 
(3.35) 
(3.36) 
(3.37) 



(an - z)x\ + a 2 ix 2 
a\ 2 xi + (an - z)x 2 
(u T + q T )xi + (v T + w T )x 2 



-(an - z)x 2 - ai 2 iEi, 
-a 21 x 2 - (an - z)xi, 



0. 



2 2 _ 7 
Ci .To — (X. 



Clearly, (I3.34|) and (|3.35ll are equivalent, and it follows from (|3.34j) and (|3.37|) that 
2 = an 



a 12 x 1 +a 21 x 2 1 

; = a n + ~ x 2 )(ai 2 xi + a 2 ix 2 ). 

xi+a; 2 a 



Next, we obtain vectors u T ,v T ,w T ,q T £ range(P) C R™ as in Case 1. Since [/' is 
an isometry, 13 . 36[) is equivalent to 



(3.38) 



(u T + q T )iri + (v T + w T )x 2 = 0. 



It follows from the analysis of cases above that rnd(^4) = 1 if and only if at 
least one of the equations (|3.33[) and (|3.38[) has a solution pair x\ , x 2 £ R satisfying 
p. lip , which proves part (i). 

If x\,x 2 satisfy (|3.33|) (resp., (|3.38p ) and (|3.1ip then y\ = x 2 , y 2 = x\ and z 
defined by (|3.14p (resp., y\ = — x 2l y 2 = —X\ and z defined by (|3.17p ) determine 
the minimal real normal completion B of the matrix A by l|3.2ip . Then the matrix 



B 



U 
1 



an 


a 12 


u 


XI 


021 


an 


V 


x 2 


T 

W 


T 

q 


S 





y\ 


!J2 





z 





0" 




A U1X1 + u 2 x 2 





1 




ViuJ + y 2 u\ z 



is a minimal real normal completion of A. Since B = B\ in Case 1 and B = B 2 in 
Case 2, this proves part (ii) of the theorem. □ 



Corollary 3.4. For a matrix A £ R" xn ; rnd(^4) = 1 if and only «/nd(^4) = 1. 
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Proof. Since we have nd(-A) < rnd(A), it suffices to prove that if nd(yl) = 1 then 
rnd(yl) = 1. Suppose that nd(A) = 1. Then, as described in Section |2~2"1 equation 
(|2.43[) has a solution x = co\(x R i, xh, xri-, x i2) € R 4 with 

(3.39) x 2 m + x 2 n > x 2 R2 + xj 2 

(see Theorem 12. 3|) . The matrix Q in this (real) case has the form 

"u T q T v T - w T 
u T + q T v T + w T 



Q 



Thus, in this case (|2.43p is equivalent to the pair of equations 

(u T - q T )x m + (v T - w T )x R2 = 0, 
(u T + q T )xn + (v T + w T )x I2 = 0. 

Since in (|3.39p either x 2 R1 > x 2 R2 or x 2 IX > x 2 2 , at least one of the equations (|3.33|) or 
(|3.38| (or equivalently, at least one of the equations (|3.9p and (|3.10p ) has a solution 
pair X\,X2 £ K with x\ > x\ (and thus, a solution pair x\,X2 £ K satisfying 
x\ — x\ = d), which by Theorem 13.31 means that rnd(A) = 1. □ 

Open problem: Is it true that for any A £ W ixn one has rnd(j4) = nd(^4)? 

As in the complex case, we will describe now a procedure (in this setting based 
on Theorem I3.3() which allows one to determine whether rnd(A) = 1 for a given 
matrix A £ M. nxn . Moreover, this procedure allows one to find all such solutions, 
and then, applying part (ii) of Theorem 13.31 all minimal normal completions of A 
when rnd(j4) = 1. 

The procedure 

Begin 

Step 1. Check the rank condition. If it holds — go to Step 2. Otherwise, stop: 
md(A) > 1. 

Step 2. Write (|3.9[) in the form (|3.33p . where u, v, w, q are defined by (|3.3ip 
and (|3.32p (see Theorem 13.31 for the definition of P, u\, and 1*2)- Find mi = 
rank(u T — q T , v T — w T ). 

Step 3. Depending on mi, consider the following cases. 

(la) mi = 0. In this case, any x\,x<z € K with x\ — x\ = d solve (|3.33p . 

(lb) toi = 1, i.e., u T — q T = ab, v T — w T = /3b, with some nonzero vector 

b e range(P) and a,f3 £ R, and additionally |a| > \(3\. In this case, ()3.33p 

is equivalent to ax\ + (3x2 — 0, and has no solutions satisfying (|3.1ip . 
(lc) toi = 1, i.e., u T — q T = ab, v T — w T = (3b, with some nonzero vector 

b £ range(P) and a,f3 6l, and additionally |a| < \(3\. In this case, (|3.33p 

is equivalent to axi + (3x2 — 0, and has the solutions 



satisfying (|3.1ip . 



(3 2 - a 2 1 * ' V (3 2 -a 2 



(Id) toi = 2. In this case, (|3.33p has no nonzero solutions, and thus has no 
solutions satisfying (|3.1ip . 
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Step 4. Write (|3. 10|) in the form (|3.38|) . where u, v, w, q are defined by 
(|3.3ip and ()3.32|) (see Theorem 13.31 for the definition of P, u\, and U2). Find 



m2 = rank(u 



T T 



+ W T ). 



Step 5. Depending on m%, consider the following cases. 

(2a) ni2 = 0. In this case, any x\, X2 € K with x\ — x\ = d solve (|3.38j) . 



(2b) m 2 = 1, i.e., u T + q T = jh, 



w = Sh, with some nonzero vector 



h € range(P) and 7, 8 € K, and additionally I7I > In this case, (|3.38|) 
is equivalent to 7x1 + 6x2 = 0, and has no solutions satisfying (|3.11[) . 



(2c) m,2 = 1, i.e. 



u 



q = jh, 



w = Sh, with some nonzero vector 



h € range(P) and 7, <5 6 K, and additionally I7I < \S\. In this case, (|3.38[) 
is equivalent to 7x1 + 5x2 — 0, and has the solutions 



Xi 



±8* 



6 2 



r 



X2 = T-f\ 



S 2 



r 



satisfying (|3.1ip . 

(2d) U12 = 2. In this case, (|3. 38[) has no nonzero solutions, and thus, has no 
solutions satisfying p. 111) . 

Step 6. rnd(A) = 1 if and only if neither of the combinations (lb)&(2b), 
(lb)&(2d), (ld)&(2b), (ld)&(2d) occur. If it does, stop: rnd(A) > 1. Otherwise, 
for each pair x\, x 2 £l obtained at Step 3, find a minimal real normal completion 
of A as described in p,13p - (|3.15p ; for each pair X\,x% G M obtained at Step 4, find 
a minimal real normal completion of A as described in (|3.15p ^ (|3.17p . 

End 

Of course, if one is interested only in checking whether rnd(A) = 1, the procedure 
can be terminated as soon as any of cases (la), (lc), (2a), (2c) occurs. 



Example 3.5. In Example [27231 



A = 



1 

1 1 

1 1 



is a matrix with real entries, and 



A A — A A 



1 

0-10 




so that the rank condition is satisfied. By Corollaries 12 . 1 1 1 and 13 . 41 rnd(A) = 1. We 
have u\ = e\, U2 — e 2 , and 



P = I — UiuJ — U2U2 — 





1 



Then, in the procedure above, u 1 



0, 



€3. Since mi 



rank(u 



T T 

q X 



0, as in Case (la), any xi,X2 € IK with x\ 



solve (|3.33[) . We have y\ = X2 and y% =x\. According to (|3.14|1 . z = 1. Therefore, 
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for any x\ £ K such that |aci| > 1, 

1 3?i 

1 1 ±y/x\ - 1 

1 1 

_±y/xf - 1 X! 1 

is a minimal real normal completion of A. We also have 

ni2 = rank(u T + q T , v T + w T ) = rank(0, 2e3) = 1, 

and as in Case (2c), h = e$, 7 = 0, 5 = 2, so that x\ = ±1 = — = = — j/i- 
According to (f3"TTT|) , 2=1. Thus, 

"10 ±1" 
110 
10 10 
=Fl 1 

is a minimal real normal completion of A. We therefore obtain the set of minimal 
real normal completions of A arising from Cases (la) and (2c). Note that the 
minimal real normal completions of A in this example are special cases of minimal 
normal completions of A as in (|2.14|) . where x\ € M: |xi| > 1, xi = ±y/ x\ — 1, and 
/i = 1, or where X\ = i, X2 = 0, and fi — ±i. 

We know from Example 12.231 that ud(^4) = 2 and that A is unitarily (and there- 
fore orthogonally) irreducible. This example shows that rnd(A) < ud(A) is possible 
for an orthogonally irreducible matrix A. 

Example 3.6. As we saw in Example l2.16l nd(-A) = ud(A) = 1, for the shift matri- 
ces A of size greater than 3. Moreover, in this case all minimal normal completions 
are actually minimal unitary completions. Since the shift matrix A has real entries, 
and 

nd(A) < rnd(A) < ud(A), 

we conclude that rnd(j4) = ud(-A) = 1, and all minimal real normal completions 
are actually minimal orthogonal completions. This corresponds to £ and p in (|2.44|) 
independently taking values 1 or — 1. 

3.3. The generic case. We will describe now the generic situation in each matrix 
dimension n. As in the complex case, a finer analysis is needed for this. However, 
in the real case our analysis is more straightforward and does not use a "heavy 
machinery" of symmetric extensions. 

For a real matrix A satisfying the rank condition, it follows from Lemma 12.71 
that the following identities hold: 

rj-\ rri rrt rr\ rr\ rr\ rr\ 

(3.40) uu = qq , vv = ww , uv = wq , uw = vq . 

Consequently, 

(u+q)(u-q) T = 0, (v+w)(v-w) T = 0, (v+w)(u-q) T = 0, (u+q)(v-w) T = 0, 

i.e., each of the vectors (u + q) T and (v + w) T is orthogonal to each of the vectors 
(u — q) T and (v — w) T . Note that the vectors u T , v T , w T , and q T belong to 
range(P) whose dimension is n — 2. 

Restricting our attention to real matrices in 9Jt„ , we now consider different values 
of n separately. 
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n = 2: In this case, vectors u T ,v T ,w T , and q T do not arise, thus necessarily 
rnd(A) = 1. This follows also from the fact the nd(A) = 1 by Corollarv l3.4l 

n = 3: In this case, vectors u T ,v T ,w T , and q T are collinear, and in view 
of (|3.40p either u T = q T and v T = w T , or u T = q T and v T = — w T . 
Thus either Case (la) or Case (2a) in the Procedure occurs. Therefore, 
necessarily rnd(A) = 1 (this follows also from Corollaries 12 . 1 1 1 and I3.4p . 

n = 4: Generically, u T ^ ±q T , and v T ^ ±w T . Since dim(range(P)) = 2, 
the vectors (u + q) T and (v + w) T are collinear and orthogonal to (u — q) T 
and (v — w) T , which are also collinear. Then both the combination of 
Case (lb) and Case (2b), and the combination of Case (lc) and Case (2c) 
(and thus, both rnd(A) = 1 and rnd(A) > 1) occur on the sets whose 
interior is nonempty in the relative topology of the manifold 9JI4. Indeed, 
the first combination occurs when we fix u T , q T and make vv T = ww T 
small enough, and the second combination occurs when we fix v T , w T and 
make uu T = qq T small enough. 

n = 5: Generically, u T 7^ ±q T , and v T ^ ±w T . Since dim(range(P)) = 3, at 
least one of the pairs of vectors (generically, only one such pair), (u + q) T 
and (v + w) T or (u — q) T and (v — w) T , is collinear. As in the case n = 4, 
for a collinear pair, both cases (b) and (c) occur on the sets whose interior 
is nonempty in the relative topology of OJI5. Thus, combinations of Case 
(lb) and Case (2d), Case (Id) and Case (2b) (where rnd(A) > 1), and 
combinations of Case (lc) and Case (2d), Case (Id) and Case (2c) (where 
rnd(A) = 1) occur on the sets whose interior is nonempty in the relative 
topology of 9Jt 5 . 

n > 6: Since dim(range(P)) > 4, the pairs (u + q) T , (v + w) T and (u — q) T , 
(v— w) T are generically linearly independent. Therefore, the combination of 
Case (Id) and Case (2d) (corresponding to rnd(A) > 1) occurs generically. 

Thus, we see that the generic situation in the real case is similar to the one in 
the complex case. 



4. Commuting completion problems 

The problem of finding commuting completions of a iV-tuple ofnxn matrices 
was raised in [3] , where a special emphasis was placed on symmetric completions of 
iV-tuples of symmetric matrices. In [8], an inverse completion (A ext , B cyA ) of a pair 
(A, B) was constructed. Namely, A cxt , B cxt by definition satisfy A ext B ext = cxl 
with a non-zero scalar a, and therefore commute. Our results from Sections [2] and 
[3] can be used to tackle commuting completion problems in the classes of Hcrmitian 
(resp., symmetric, or symmetric/ antisymmetric) pairs of matrices. 



4.1. The commuting Hermitian completion problem. Let (A\, A2) be a pair 
of Hermitian matrices of size n x n. Wc define the commuting Hermitian defect of 
A\ and A2, denoted chd(Ai, A2), as the minimal possible p such that there exist 



commuting Hermitian matrices B\ = 



Ai 



* 



and B2 = 



A, 



of size (n+p) x 
chd(A 1 ,A 2 )) 



(n + p). We call such a pair (B\, B2) of size (n + chd(Ai, A 2 )) x (n 
a minimal commuting Hermitian completion of (A\, A2). 

Since (B%, B2) is a commuting Hermitian completion of a pair (A\, A 2 ) of Her- 
mitian matrices if and only if B = B\ + iB 2 is a normal completion of A = A1+1A2 
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and therefore chd(Ai, A 2 ) — nd(Ai + i4 2 ), the results from Section l2~2l allow one 
to check whether chd(Ai, A2) = 1, and when this is the case — to construct all 
minimal commuting Hermitian completions of (A\, j4 2 )- For example, Theorem l2.3l 
yields the following. 

Theorem 4.1. Let A±,A 2 G C nx ™ be Hermitian. 

(i) chd(j4i, A2) = 1 if and only £/ rank(j4ij4. 2 — A2A1) = 2 and the equation 

(4.1) PA 1 {t 1 u 1 -t 1 U2) = iPA 2 {t2U 1 +t2U 2 ) 
has a solution pair ti,t2 G C satisfying 

(4.2) Re(?ii 2 ) = d. 

Here Ui,u 2 G C™ are the unit eigenvectors of the matrix 2i(A\A2 — A2A1) cor- 
responding to its nonzero eigenvalues Ai = d(> 0) and A2 = —d, and P = 
In — u\u\ — U2U2 is the orthogonal projection of C™ onto null(j4iA 2 — A2A1). 

(ii) 7/chd(vli, A2) = 1, t\ andt2 satisfy (|4.1|) and (14. 2|) . and /1 G T is arbitrary 
then the pair (B\,B2) of matrices 

A x §(t 2 ui +I2U2)' 

J2 {tml + t 2 U2) Zi 



(4.3) 



Bi = 



(4.4) 



(4.5) 



B 2 = 



A, 



1 - hu 2 ) 

Z2 



is a minimal commuting Hermitian completion of (A\, A2) 
1 



Here 



u\Aiui — — (lm(t 2 U2 J ^2Ui) + Re(^ii2M2^i u i)) 



(4.6) 



z 2 = u\A 2 ui 



- (lm.{t\u*2AxUi) - Re(iii 2 u2^2Ui)) • 



All minimal commuting Hermitian completions 0/(^1,^2) arise in this way. 

Proof. Letting A = A x + iA 2 , we observe that A* A - AA* = 2i(A x A 2 - A%A\). 
It is straightforward to verify that, under the change of variables ti — x\ — X2, 
t2 = X1+X2, condition (|2.1ip in Theorem 12.31 is equivalent to condition (|4.1[) . 
(|2.12p is equivalent to |42)l . B x and B 2 defined by l|4.3p and (j4~4)) are Hermitian 
and such that B = B\ + i£>2 is as in (|2.14p , z\ and Z2 defined by (|4.5|) and ()4.6j) 
are real and such that z = z\ + iz2 is as in (|2.15|) . Thus, this theorem is equivalent 
to Theorem [O □ 



4.2. The commuting completion problem in the class of pairs of symmet- 
ric and antisymmetric matrices. Let A\ = G R nxn and A 2 = -A\ G R nxn . 
It is natural to ask what is the minimal p such that there exist commuting matrices 



B x =Bl = 



G R(™+p) x (™+p) and B 2 = -BT = 



Ao 



o(n+p)x(n+p) 



Such a pair (B\,B2) is a minimal commuting completion of (4 1; 4 2 ) in the class 
of pairs of symmetric and antisymmetric matrices. Since (Si,^) is a commuting 
completion of (At, A 2 ) in this class if and only if B = B\ + B2 is a real normal 
completion of A = At + A2, our results from Section [3] can be restated in terms of 
pairs of matrices in this class. We omit the details, since the reasoning is similar 
to the one in Section |4~T1 



MATRICES WITH NORMAL DEFECT ONE 



31 



4.3. The commuting symmetric completion problem. In this section, we 
consider the commuting completion problem in the class of pairs of symmetric 
matrices. This is a special case of the problem raised in Degani et al. [2j (see the 
first paragraph of Section [4j for N = 2. The authors of [3] presented an approach 
to n-dimensional cubature formulae where the cubature nodes are obtained by 
means of commuting completions of certain matrix tuples. While their commuting 
completion problem is stated in a certain subclass of tuples of symmetric matrices, 
some observations were also made for the problem in the whole class. In particular, 
the question on the minimal possible size of completed matrices was accentuated 
as important. 

Let Ai = A[ € K" x ™ and A 2 = A$ £ W lXn . We define the commuting sym- 
metric defect of A\ and A 2 , denoted csd(Ai, A 2 ), as the minimal possible p such 

\A\ * 

that there exist commuting symmetric matrices B\ = 



,B 2 = 



A 



R (n+p)x(n+ P )_ Such a pair (B U B 2 ) of size (n + csd(A 1 ,A 2 )) x (n + csd(Ax, A 2 )) is 
a minimal commuting symmetric completion of the pair {A\,A 2 ). 

We note that (B%, B 2 ) is a commuting symmetric completion of a pair (Ax, A 2 ) of 
real symmetric matrices if and only if B = B\ + iB 2 is a normal, and simultaneously 
complex symmetric, completion of A = A\ + iA 2 . We also observe that a priori 

(4.7) csd(Ai, A 2 ) > chd(Ai,A 2 ). 

Open problem. Is it true that for any pair (A\,A 2 ) of real symmetric matrices 
one has cs&(A\,A 2 ) = ch.d(Ax,A 2 )7 

This problem is somewhat similar to the open problem stated in Section [3~2l The 
latter actually asks whether a minimal normal completion of a real matrix can be 
chosen to be real, while the former is equivalent to the question whether a minimal 
normal completion of a complex symmetric matrix can be chosen to be complex 
symmetric. 

The following theorem shows that, for a pair (Ax, A 2 ) of real symmetric matrices, 

csd(A x ,A 2 ) = 1 <=^ chd(^i,A 2 ) = 1, 

which motivates the open problem above. Moreover, this theorem actually shows 
that if csd(Ai, A 2 ) = 1 then the set of all minimal commuting symmetric comple- 
tions (Bi,B 2 ) of (Ai,A 2 ) can be obtained by putting in Theorem 14. II u-> = u\ and 
/i = l. 

Theorem 4.2. Let Ai,A 2 S M. nxn be symmetric. 

(i) csd(Ai, A 2 ) = 1 if and only if ia,'nk(AiA 2 — A 2 A\) = 2 and the equation 

(4.8) PA X Im(tiUi) = PA 2 Re(t 2Ul ) 
has a solution pair tx,t 2 G C satisfying 

(4.9) Re(iii 2 ) = d. 

Here u\ € C" is the unit eigenvector of the matrix 2i(AxA 2 — A 2 Ax) corre- 
sponding to its eigenvalue X\ = d(> 0), and P = I n — uxu\ — Hxuf . 

(ii) 7/csd(Ai, A 2 ) = I, tx and t 2 satisfy (14. 8|) and (|4.9p then the pair (B 1} B 2 ) 
of matrices 

Ax Re(t 2 ux) 
Re(t 2 ux) T zx 



(4.10) Bx = 
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(4.11) 



A 2 Im(tiMi) 
[lm(^iUi) T z 2 

is a minimal commuting symmetric completion of (A\, A 2 ). Here 



(4.12) 



(lm(t|uf A 2 ui) + Re(tii 2 w 



ZAxUi)) 



and 



(4.13) 



z 2 = u\A 2 ui 



(lm(^uf Aim) - Re(tit 2 u^ A 2 ui)) . 



All minimal commuting symmetric completions of {A\, A 2 ) arise in this way. 

Proof, (i) By (077]), if csd{A u A 2 ) = 1 then chA(A x ,A 2 ) = 1. Therefore, by 
Theorem 14. 1[ r&nk(AiA 2 — A 2 Ai) = 2 and equation (|4.1[) has a solution pair 
ti,t 2 £ C satisfying (|4.2p . If u\ is the unit eigenvector of the Hermitian matrix 
2i{A\A 2 — A 2 A{) corresponding to its eigenvalue Ai = d(> 0) then u\ is the unit 
eigenvector corresponding to the eigenvalue A2 = —d. Thus, we can choose in The- 
orem |4~T1 it? = Ttx- Then P = I n — u\u\ — Tiiu[ is a real n x n matrix, and equation 
dHH) becomes 

Conversely, if rank(AiA2 — A 2 Ai) = 2 and equation (|4. 8|) (which is equivalent 
to (|4.ip in our case) has a solution pair t\,t 2 £ C satisfying (|4.9p (= l|4.2p ) then 
by Theorem 14. li chd(Ai,^42) = 1- For any such t\,t 2 the corresponding minimal 
commuting Hermitian completions (P>i,P> 2 ) of (Ai,A 2 ) have the form (|4.3| - (|4.4I1 . 
We observe that since u 2 =u\, the matrices P>\ and B 2 are real symmetric if and 
only if /j = 1 or fi = —1. Consequently, csd^i,^) = 1. 

(ii) If ti,t 2 £ C satisfy g^])-gl]) then so do t[ = -t x and t' 2 = -t 2 . There- 
fore, we do not miss any minimal commuting symmetric completions of (Ai,A 2 ) 
if in (|4.3p ~ (|4.4p we choose t x ,t 2 as above and fix /i = 1. Finally, since (Bi,B 2 ) 
constructed in Theorem 14. II with /1 = 1 has in our case the form (|4.10p - (l4.1ip . and 
(H3D-(|MD become (|4TT2]) - (|4~r3l) . this completes the proof. □ 

Remark 4.3. The procedure for checking whether csd(Ai,^42) = 1, and if this is 
the case — for finding all minimal commuting symmetric completions of a pair of 
symmetric matrices (Ai, A 2 ), can be obtained as the specialization of the procedure 
mentioned in Section 14.11 (which, in turn, is based on the procedure from Section 
12. 2p by setting u 2 = u\ and [i — 1 (see Theorem 14.21 and its preceding paragraph) . 



In the 1980s the use of quantum systems as computing devices started to being 
explored. The idea gained momentum when Peter Shor [JO] presented a quantum 
algorithm for factoring a large composite integer N that was polynomial in the 
number of digits in N. An excellent overview article on the subject of quantum 
computing is pQ. 

The separability problem occurs when a quantum system is divided into parts. 
For convenience we consider a bipartite system. The state of the system is described 
by a density matrix M, a positive semidefinite matrix with trace 1. A state is called 
separable when it can be written as a convex combination of so-called pure separable 
states, i.e., p = Yli=iPi V'iV'i* ® ^i^X where ipi and (pi are (nonzero) state vectors in 
the spaces corresponding to two parts of the system, and pi > 0. When £ C m 



5. The separability problem 
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and <j>i £ C", the matrix p is called m x n separable. The number k is referred to 
as the number of states in the representation. 

The problem whether a given state is separable or entangled (= not separable) 
may be stated as a semi-algebraic one, and is therefore decidable by the Tarski- 
Seidenberg decision procedure [2] . As it turns out though, the separability problem 
scales very poorly with the number of variables and these techniques are in general 
not practical. In fact, the separability problem in its full generality has been shown 
to be NP-complete [5]. 

As a consequence of the results of Section [2] we can state a new result for the 
2 x n case. Thus we are concerned with matrices 



(5.1) 



M 



A B* 
B C 



> 0. 



Notice that if M = Ylt=i p * ® Qi with p i e 
C nx " positive semidefinite, then M = Y2i=x P T 
well. One easily sees that 



2x2 positive semidefinite and Qi € 
Qi is positive semidefinite as 



(5.2) 



M 



A 

B* 



> 0. 



Thus for (|5.1[) to have a chance to be 2 x n separable we need (|5.2I) to hold. This 
is referred to as the "Peres test"; see [3]. As was observed by several authors, the 
2xn separability problem for (|5.ip can easily be reduced to the case when A = I; 
see, for instance, Proposition 3.1 in [12] . Using Theorem 3.2 in [12], which connects 
the separability problem to the normal completion problem, we can now state a 
method for checking separability of (|5.ip in the case when rank(M) = rank(M) = 
rank(A) + 1. 



Theorem 5.1. Let B, C e 



be such that 



(5.3) 

and suppose that 

(5.4) 

Write 



M 



In B* 
B C 



> 0, M 



In B 

B* C 



>0, 



rank(M) = rank(M) = 



C 



BB*=xx*, C-B*B = yy* 

for some vectors x, y G C". Then M is 2 x n separable if and only if x 7 y, B*x, By 
are linearly dependent. In this case, the minimal number of states in a separable 
representation of M is n 



1. 



Proof. First notice that B*B — BB* = xx* — yy* . 

Suppose that x, y, B*x, By are linearly dependent. Then by Theorem 
exists a normal matrix 

B vx 



there 



N 



y 



where |^| = 1. But as (vx)(vx)* = C — BB* it follows from Theorem 3.2 in [12] 
that M is 2 x n separable, and that the minimal number of states in a separable 
representation of M is n + 1. 
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Conversely, suppose that M is 2 x n separable. By Theorem 3.2 in [12] there 
exists a normal matrix 

B S~ 
T P 



N = 



so that BB* + SS* < C. But then SS* < xx* and thus S = xv* with \\v\\ < 1. 
Also B*B + T*T = BB* + SS* < C, and thus T*T < yy* yielding T = yw* with 
||w|| < 1. In addition, BT* + SP* = B*S + T*P. In particular, range(BT* + 
SP*) = tsmge{B*S + T*P). Note that range(ST* + SP*) C span(By,x) and 
range(B*S' + r*P) C span(i?*a;, y). But then it follows easily that x, y, B*x, By are 
linearly dependent. Indeed, if BT* + SP* = B*S + T*P ^ 0, then span(B?/, x) and 
span(B*cc, y) must have a nontrivial intersection, and if BT* + SP* = B*S+T*P = 
0, then spa,n(By, x) and span(B*x,y) are both at most one dimensional. □ 

We can now provide a new proof of the following result by Woronowicz |13) . 

Theorem 5.2. Let A, B,C be n x n matrices with n < 3, so that 



(5.5) 



M 



A B* 
B C 



> 0, M 



A 

B* 



> 0. 



Then M is 2 x n separable. 

We will use a result by Hildebrandt which we quote without proof. 

Lemma 5.3. [6J Lemma 2.6] Let K be a convex cone in a real vector space TL of 
finite dimension N , and let C QTi. be a subspace of dimension n. Let K' = K D C 
and x the generator of an extreme ray in K' . Then the minimal face in K containing 
x has dimension at most N — n + 1 . 

Notice that if we consider the cone PSD n of n x n complex positive semidefinite 
matrices, then the minimal face containing M > 0, is the cone F = {GCG* : C £ 
PSD k } 7 where M = GG* with null(G) = {0}, and k = rank(M). In particular, 
the real dimension of this minimal face is (rank(M)) 2 . 



Proof of Theorem \5.SX Since the case n < 3 can be embedded into the case n = 3, 
we will focus on the latter. As the 2 x n separable matrices form a convex cone, 
it suffices to prove the result for pairs (M, M) that generate extreme rays in the 
cone of pairs of matrices as in (|5.5p . If we apply Lemma 15.31 with the choices of 
K = PSDq x PSDq and C the subspace 



B* 
C 



A 

B* 



in the (real) vector space of pairs of Hcrmitian matrices of size 6x6, then K' 
is the cone of pairs of matrices as in (|5.5p . By Lemma 15.31 the minimal faces 
in K containing extreme rays of K' cannot have dimension greater than 72 — 
36 + 1 = 37. However, the minimal face in K containing (M, M) (which gener- 
ates an extreme ray in K') has dimension (rank(M)) 2 + (rank(M)) 2 , and hence 
the vector (rank(M), rank(M)) £ R 2 lies in the closed disk of radius V37 cen- 
tered at the origin. This now gives that either min{rank(M ), rank(M )} < 3 or 
max{rank(M), rank(M)} = 4. Next, as in Proposition 3.1 in 12] we can assume 
that A = I. If now min{rank(M), rank(M)} < 3 we have that C = BB* = B*B, 
and thus B is normal, which yields by Theorem 3.2 in [T2] that M is 2 x n separable. 
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On the other hand, if max{rank(M), rank(M)} = 4 we can conclude by Theorem 
15.11 that M is 2 x n separable (as 4 vectors in C™ are always linearly dependent 
when n < 3). □ 

It should be noted that the original statement of Woronowicz is formulated in 
the dual form: if $ : C 2x2 — * C" xn is a positive linear map (thus $(PSD 2 ) C 
PSD n ) and n < 3, then $ must be decomposable. That is, $ must be of the form 
= Eli RiMRf + ELi S t M T S*. 

Acknowledgment. The authors wish to thank Dr. Roland Hildebrand for pointing 
out how Theorem 15. ll in conjunction with his results, leads to a proof of the result 
by Woronowicz. 
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